Multichannel field-effect spin barrier selector 
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We have studied spin carrier dynamics under full spin-orbit coupling. The anisotropy of disper- 
sions for independent circular spinor polarizations is explored as a possible vertical multichannel 
voltage controlled spin-filter. Small voltage variations are found to select the current polarizations 
in a resonant tunneling geometry. 
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The advent of spintronics has resulted in the study 
and design of spin manipulated devices used as infor- 
mation processors, quantum computing elements, spin- 
polarized diodes, spin-valve read headsji and electro- 
optical modulators^ to name just a few. The spin-orbit 
(SO) interaction, due to spatial asymmetry in zincblende 
lattices or sample design, doping profile or applied gate 
voltages, plays a fundamental role on these systems? 3 *^ 
especially on spin-dependent tunnelling. The SO interac- 
tion couples the electronic momentum to the spin degrees 
of freedom, and lifts the spin degeneracy for structures 
fabricated in zincblende materials. The narrower the en- 
ergy gap of the host material, the stronger these effects 
will appear on the transport and optical properties^ 

We present here a detailed analysis of the complexity 
of bulk electronic structure, derived from kinetic energy 
plus the full SO Hamiltonian including Bychkov-Rashba 7 
and all three Dresselhaus^ contributions. Then, we study 
SO effects on spin-polarized current in double-barrier res- 
onant (DBR) devices, and analyze how the anisotropy on 
the spin orientation and polarization of spinor states can 
be used to select and optimize channels for vertical trans- 
port in the system. 

Recent publications have proposed the manipulation 
of linear Rashba and Dresselhaus terms in the SO 
interaction^ The cancellation of linear terms leads to 
drift-diffusive lateral transistors, in contrast to the bal- 
listic operation of the Datta-Das device^ Perel et aU& 
have also reported on the evolution of the spin orienta- 
tion with 2D transverse momentum, and explored the 
tunnelling through a single barrier with different III-V 
materials, driven only by the linear k\\ Dresselhaus term. 
A time-dependent spin manipulation scheme has been 
recently proposediii 

Let a = + (— ) label the spin-up (spin-down) state. 
The carrier dynamics, driven by the kinetic energy plus 
full SO Hamiltonian, is dictated by 

e *u-P. 
(1) 

In this equation, the diagonal terms are H aa — 
- (h 2 /2m*)d 2 /dz 2 + (h 2 k 2 /2m*) - iaa 2 (fe|| , <p)d/dz, while 
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the off-diagonal term is H + _ = /3n(kn , ip) + 03 (fcii , p) + 
ai(k\\,<p)d 2 /dz 2 , with H |_ = H\__. In this con- 
text, k|| = (k x ,k y ) = (fc||,<^) is the in-plane wave 
vector, p — (x, y) is the electron position in the xy 
plane, cki (k\\ , p) = aok^e lip , a>2(k\\,(p) — aok 2 cos2ip and 
ct3(k\\,<p) = a k^(e lip — e~ l3ip )/4 are, respectively, the 
linear, quadratic and cubic Dresselhaus SO terms. Fi- 
nally, P R (k h p) = (3 (-dV R /dz)k\\e l ^i 2 -^ is the linear 
Bychkov-Rashba contribution under any potential profile 
Vr(z), whereas F±(z) are the spinor components. Here, 
«o and /?o are the SO material constants. 12 

First let us discuss bulk dispersions where k z is a 
good quantum number for the —id/dz operator. Due 
to the strong coupling between spin s and linear mo- 
mentum k = (kii,fe z ) degrees of freedom, the Dressel- 
haus (or bulk inversion asymmetry, BIA) and Bychkov- 
Rashba (or surface inversion asymmetry, SIA) Hamil- 
tonians create special sectors and anisotropics on the 
spinor phase-space. For given k|| and applied electric 
field F = —(l/e)dVn/dz, the bulk eigenvalues of the full 
Hamiltonian in Eq. ^ are 



£±0*1,*.) = — (fcrf + K) ± Ji(k h k z ) + *(k,,), (2) 
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where j(k.u,k z ) 



a^knk 2 cos(4p) 



* n,||i 



2a pk 2 k 2 sm(2(p) and <5(k||) = /3 2 k 2 + a 2 k\ sin 2 (2<p)/4 + 
ao/3fc| sin(2</?). Here ip is measured from the (1,0) 
crystalline axis on the xy plane and (3 — (3o(—dVn/dz). 

Notice that from the four degeneracies, two for the 
momentum (±fc z ) and two for spin (±), one can con- 
struct only two independent spinor states forming the 
Kramers doublet. Without SO, there are two degener- 
ate states, well known as linear spin polarized, where the 
spin-up and spin-down spinors display isotropic parabolic 
dispersions for both k\\ and k z (first term in Eq. [5J. 
However, with SO there are two orthogonal Hilbcrt sub- 
spaces, the circular spin polarized states, with highly 
anisotropic dispersions when the orbital angle, pi, changes 
clockwise (cr + ) and counterclockwise (er - ), respectively. 
The gap measuring the spin-splitting is E a g (k^,p,k z ) — 
[E + (fcy , p, k z ) — E-[k\\, p, k z )], or twice the square root 
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FIG. 1: Angular anisotropies and bulk dispersions for de- 
polarization at special angles within the interval < tp < 
180°. Notice spin polarization in different sectors of k z . Pan- 
els at tp — tp a c = 34.842° and at tp = 145.158° show no spin 
gap. 



in Eq. |5| At k z = 0, this yields a critical angle tp a c = 
-0.5sin- 1 [-2/3/(a /cj|)] where E°(k\\ , tp a c , k z ) = 0. This 
result depends on the field strength and on the ratio be- 
tween SO constants, and defines a maximum value for 
given by fen = y/— 2/3/ao- There are two symmetry 
regimes according to which SO mechanism is dominant: 
(a) SIA regime, or Rashba dominant, for k\\ < fefj , and (b) 
BIA regime, or Dresselhaus dominant, for k\\ > k?i. These 
properties are valid for any zincblende material and here 
we will fix our attention on InSb, due to its large SO 
parameters^ These states and their properties are rem- 
iniscent of circularly-polarized light travelling through a 
crystal, as we shall discuss later. It is convenient to de- 
fine a new zero for the angle at the (1,-1) crystalline 
direction in the xy plane, corresponding to a rotation 
tp => ip — 45°. We adopt this from now on. 

Figure ^ shows different cuts on the energy surface, 
for the BIA regime and o~ + polarization, calculated for a 
fixed value of fen = 3, and different values of tp when 
a uniform field, Fq = 5, is applied. For InSb— and 
the given Fo, the critical values are tp° = 34.842° 
and fcjj = 1.768. The panel for tp = shows the non- 
parabolic spin-splitted dispersions for spin-up and spin- 
down branches. For a given energy inside the gap re- 
gion, the spin-up states have imaginary values for k z 
along a real line, shown as dotted green lines in the 
panel with tp = 90°, and these states do not propa- 
gate. The next panel with tp — tp a c shows the situation 
where E g {k\\,tp° ,k z = 0) = 0. Just above this angle, 
at tp = 35°, the gap opens again and an inverted group- 
velocity layer of spin-up states (inside dotted blue-box) 
appears between critical fc z -values, ±fc c ij where the z- 
group-velocity vanishes. The existence of this branch is 
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FIG. 2: Spin-splitting anisotropy for SIA regime (at feii = 
1.768 = fefj, panel a) and BIA regime (at fcy = 3 > fcjj, panel 
c) produced by full SO (black/red lines) and BIA terms only 
(green lines) in bulk InSb. Panels b and d: Twofold and four- 
fold symmetries for spin-splitting gap in the SIA (b) and BIA 
(d) regimes, at different values of feii . Panel e: Comparison 
between Rashba and Dresselhaus interactions, for panel c con- 
ditions; red (green) lines for BIA+SIA (BIA) terms. Panel /: 
Spin-splitting anisotropy of a + energy surfaces, away from 
zone-center. 



determined partly by the quadratic 0:2 (fell, <f) BIA term 
in Eq. ^ For k z in the vicinity of these critical points we 
have -E+(k|| , k z ) ~ E +S (k z ±k c i) 2 , with E and con- 
stants. Inside the dotted blue box the state has inverted 
group-velocity and negative effective mass, such that a 
carrier will travel to =Fz-axis for ±fc 2 values of the propa- 
gating wave vector. In the higher spin-up and spin-down 
branches outside ±fc c i, the carriers travel with normal 
group- velocity. By increasing the angle to tp — 45°, the 
double-valley region disappears, with an increasing gap 
that reaches its maximum value at tp = 90°. The panels 
with 90° < tp < 180° just complete half of the periodicity 
of the cr + Hubert subspace. 

The dispersions for o~ polarization, in this BIA 
regime, have identical features as in Fig. ^ except that 
they require exchange between spin-polarizations in the 
sectors. The Hilbert subspaces form the degenerate 
Kramers doublet for circular polarizations. These a~ 
and o~ + spinors are characterized by their behavior un- 
der the time-reversal operator for zincblende symmetry, 

T = —iijyZl , where —icr y flips spin and the complex con- 
jugation (inversion) C (i) flips momentum (position). 

Figure shows the angular gap variation, for both 
circular polarizations, at the critical situation where 
fc|| = fcjj — 1.768, for Fo = 5. Black and red lines are 
the combined results for BIA+SIA terms whereas green 
lines show the effect of BIA terms acting alone. It be- 
comes clear that BIA+SIA terms act in-phase and out- 
of-phase in each sector. This alternating change in phase 
between the SO Hamiltonian terms leads to a Rashba 
regime with Ciy symmetry, as shown in Fig. |2t> for dif- 
ferent fc|| < fcjj. In this regime, the Hamiltonian in Eq. ^ 



3 




(10° cm"' ) 



FIG. 3: 2D-subband dispersions, solutions of Eq.0 in isolated 
InSb QW with L w — 200 A and uniform electric field Fo = 
100. Different anisotropics induced by the full SO terms result 
in varying in-plane masses, m*(k\\,(p), and spin-splittings for 
each subband. Red (black) lines for down (up) spin. 



produces a smooth change from s- to p-like symmetry on 
the electronic structure at the zone-center as ku increases. 



Figure yj; shows the angular gap variation in the Dres- 
selhauss regime (black and red lines for fc|| = 3 > A;Jj, for 
Fo = 5). The effect of BIA terms alone is shown in green 
lines. Here also the full SO Hamiltonian acts in-phase 
and the overall symmetry is changed from a twofold C2V 
to a fourfold C±y. As shown in Fig. 01, this modification 
induces smooth changes on the electronic structure at the 
zone-center, from p- to d-like symmetry as ku increases. 
This is also at the origin of the secondary gap and the 
double-valley sector near the critical angle. The result- 
ing symmetries show that it is essential to consider all the 
three BIA terms in the Hamiltonian, so that their cumu- 
lative effects yield the correct electronic structure and 
eigenstates. Observe that at given electric field there are 
only specific angles where the Rashba and Dresselhaus 
effects cancel each other. 

Figure 12k compares the combined effects of BIA and 
SIA terms (red lines) on the spin-splitting, at the zone- 
center, for the same conditions shown in Fig. |2t and the 
red line in Fig. (notice scale changing) . It is clear that 
their effects are added along the (±1, ±1) directions of 
maximum spin-splitting and the BIA terms (green lines) 
dominate along (±l,=pl) directions where the secondary 
gap appears. 

If we move away from the zone-center, the complexity 
increases. In Fig. U, we start with the critical situation 
displayed in Figs. |2t an d 03- The cuts at increasing 
values of k z show new gap maxima along the (±1, 0) and 
(0, ±1) directions. These curves are a reflection of the 
complex geometry of the 2D Fermi surface at different 
electron concentrations. 

Evidence of strong anisotropy on FIR light propagat- 
ing through n-type InSb bulk samples was studied by 
Gopalan et ali& They have shown why the magneto- 
transmission experiments of Dobrowoslka et al^ in the 
Voigt and Faraday configurations, displayed maximum 
spin-flip and maximum transmissivity along the (±1, ±1) 
and (±1, =Fl) directions. The electronic structure we re- 
port here arises from the same complex lattice symmetry. 

Having shown the complexity of the bulk spinor en- 
vironment under the full SO Hamiltonian, let us now 



explore these anisotropies via the vertical transport of 
carriers in a DBR structure with interfaces [thickness] at 
positions zi (£ = 1,...,4) [di — (2^+1 — zt)\ and profile 
Vq(z). The layer t = (l = 4) is the emitter (collector) 
contact. We use the scattering matrix technique to cal- 
culate the transmissivity and reflectivity for the system. 
For a given incident energy, E, we want to construct the 
incoming and outgoing spinor states that tunnel through 
the 2D states in the layer 1 = 2. In Fig. [3] we show 2D 
dispersions derived from Eq. ^ in an isolated quantum 
well under applied electric field. We observe that the 3D 
anisotropies are transferred to the subband dispersions, 
through anisotropic in- plane effective masses m*(k\\,tp), 
as well as spin-splittings of each subband. 

We construct each component of a spinor state, 
F(p, k|| , z) in Eq. as a linear combination of bulk 

solutions 15 = Ea(+^)[ a f( k \\) F A+kza)e +tk ^ z ~^ 

+ bfik^FA-k^e-^-^l where of(k||) (6f(k|,)) 
is the amplitude of incoming (outgoing) spin-polarized 
waves at any given interface zi, and ±fc ZCT are real roots 
of E a (k\\,k z ) — E, for the carriers travelling along the 
z-axis. The dependence on k|| defines the open channels 
in phase-space. The current operator from Eq. ^ is 



Jz(k\\,k z ) 



J++ J+- 



(3) 



where 3 aa = hk z /m* + aa2(ku,(p)/h represents 

the spin-conserving polarized current, and J_| = 

—2m*aii(Jc\\,(p)k z /h the spin-flip mechanism. The ver- 
tical transport properties are also determined by the 
quadratic (diagonal) and linear (spin-flip) BIA SO 
terms. The boundary conditions require that both spinor 
F(p, kii, z = z{) and flux J z F(p, k|| ,z — z{) be continuous 
across each interface, z\. Resonant tunnelling for a given 
circular polarization implies that the matching conditions 
must be valid for each k|| . 

The elements of the scattering matrix S are the spin- 
conserving (ti and tZ) and spin- flip (t_ and fT) trans- 
mission, and the equivalent reflection coefficients r±, as- 
sociated to the incoming and outgoing wave amplitudes 
t°, = a% /oq, for detected (a 1 ) and incident (a) polariza- 
tions. For any applied potential V, the partial transmis- 
sivity can be calculated^ from the ratio between detected 
J°"' and incident J" 1 currents as T a ^, a * (E, V, k|| ) = 
Ke[(t°,)n°,{F^J+ a: \F+,) z=z J (F+\J+ a \F+) z=Z0 \. With 
the presence of spin-flip processes, the partial reflec- 
tivity is calculated from the general flux conservation, 
EAT*^*>{E,V) + R a ^ a ,(E,V)] = 1, for each inci- 
dent polarization. We consider vertical transport along 
the z-direction (001) of a symmetric InSb-InAlSbiS DBR 
[l b = 100 A (barrier) and L w = 200 A (well)]. 

Figure 0] shows calculated partial transmissivities in 
the <r + polarization. In principle, at a given incident en- 
ergy E, there will be transmissivity peaks at each spin- 
polarized QW level. As mentioned above, spin-up states 
cannot propagate when E is in the gap-region since the 
k za roots of E a {\n\\,k z ) = E are imaginary, as shown 
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FIG. 4: Transmissivity for a spinor in InSb-InAlSb dou- 
ble barrier with lb = 100 A and L w = 200 A, calculated 
for fen = 1, Fo = 30 and ip = 35° (90°) for an incident en- 
ergy above (inside) the gap, shown in blue (green) lines. The 
up/down arrows on left (right) indicate incidence (detection) 
polarizations in the DBR. Curves are shifted upwards for clar- 
ity, and small differences in peak positions of the two highest 
curves are due to distinct values of ip. Notice different ampli- 
tudes of each resonant peak at spin-polarized 2D levels inside 
QW. For energies inside the gap, only one spin propagates. 

in panel ip = 90° of Fig. ^ F° r E outside the gap- 
region all polarizations propagate. The four transmis- 
sivity curves in Fig. 0] for an energy above the gap (blue 
lines) show double peaks at the quantum well resonances, 
for <7 + -polarization. These same features occur for (re- 
polarization at equivalent angles. However, when the in- 
cident energy is inside the gap (green lines), notice that 
only the spin-down incident polarization can propagate, 
while the spin-up channels result in vanishing transmis- 
sivity. The resulting current arriving at the collector 



will have Stark shifted double peaks in voltage for each 
quadrant^ Since a + - and er _ -polarizations belong to or- 
thogonal Hilbert subspaces, the vertical current produced 
by these states (for energies above the gap) define eight 
independent emitter-collector transmission channels that 
may be explored as voltage controlled spin-filters. In fact, 
polarized currents with desired polarization (±) can be 
produced with small voltage variations. Notice also that 
for energies inside the gap the total current is, de facto, 
spin-down polarized. 

In conclusion we have studied how anisotropics of 
spinor states and dispersions of the full SO Hamiltonian 
in orthogonal Hilbert spaces, with opposite circular spin- 
polarized configurations, could enhance the formation of 
spin- filters. The possibility to tune emitter Fermi level 
and collector detectors requires designed magnetic con- 
tact masks on the DBR diode, a technology in full devel- 
opment nowadays^ This would help exhibit and explore 
the complexity of the eight channels produced by the full 
SO Hamiltonian. We have shown that this new struc- 
ture can work as multichannel voltage controlled spin- 
filters, similar to semimagnetic materials^ We antici- 
pate that new designed experiments will soon be able 
to probe the gap structure and anisotropies we predict. 
Finally, the possible scattering mechanisms within these 
two orthogonal spinor configurations must be analyzed 
under the unusual situation where spin and momentum 
coordinates are strongly coupled. Therefore, new car- 
rier spin-flip relaxation times, tsf, proportional to the 
reciprocal of spin-splitting energy^ Eg (k\\, (p c , k z ), and 
momentum scattering times, t p , can be measured since 

Tgp = t p E " ( fc ii^ c,fc *) 1 2, The very large values for 

Eg(k\\, tp c , k z ) (~ 5 — 10 meV) can produce devices with 
tsf many orders faster than t p . 

Work is supported by FAPESP, CNPq, CAPES, US 
DOE grant no. DE-FG02-91ER45334, and the CMSS 
Program at OU. 



1 Semiconductor Spintronics and Quantum Computation, 
edited by D.D. Awschalom, D. Loss, and N. Samarth 
(Springer, Berlin), 2002, and references therein. 

2 S. Datta, and B. Das, Appl. Phys. Lett. 56, 665 (1990). 

3 C.F. Destefani, Sergio E. Ulloa, and G.E. Marques, Phys. 
Rev. B 69, 125302 (2004). 

4 H.C. Liu, J. Li, G.C. Aers, C.R. Leavens, M. Buchanan, 
and Z.R. Wasilewski, Phys. Rev. B 51, 5116 (1995). 

5 J. Nitta, T. Akazai, and H. Takayanagi, Phys. Rev. Lett. 
78, 1335 (1997). 

6 G.E. Marques, and L.J. Sham, Surf. Sci. 113, 131 (1982). 

7 Y.A. Bychkov, and E.I. Rashba, J. Phys. C: Solid State 
Phys. 17, 6039 (1984). 

8 G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

9 X. Cartoixa, D.Z.-Y. Ting, and Y.-C. Chang, Appl. Phys. 
Lett. 83, 1462 (2003); J. Schliemann, J.C. Egues, and D. 
Loss, Phys. Rev. Lett. 90, 146801 (2003). 

V.I. Perel, S.A. Tarasenko, I.N. Yassievich, S.D. Ganichev, 



V.V. Bel'kov, and W. Prettl, Phys. Rev. B 67, 201304(R) 
(2003). 

11 E. I. Rashba, and Al. L. Efros, Phys. Rev. Lett. 91, 126405 
(2003). 

12 Parameters for InSb/Ino.72Alo.2sSb: ao = 500 A 2 ; /3o = 
160 eVA 3 ; m* = 0.0143m /0.043m ; offset: Vo = 300 
meV. Units used in the paper: meV for energy, 10 6 cm -1 
for wave vectors, kV/cm for electric field. 

13 S. Gopalan, J.F. Furdyna, and S. Rodriguez, Phys. Rev. B 
32, 903 (1985). 

14 M. Dobrowolska, Y.-F. Chen, J.F. Furdyna, and S. Ro- 
driguez, Phys. Rev. Lett. 51, 134 (1983). 

15 A.C.R. Bittencourt, A.M. Cohen, and G.E. Marques, Phys. 
Rev. B 57, 4525 (1998). 

16 E.I. Rashba, Phys. Rev. B 62, R16267 (2000). 

17 A. Slobodskyy, C. Gould, T. Slobodskyy, C.R. Beker, 
G. Schmidt, and L.W. Molenkamp, Phys. Rev. Lett. 90, 
246601 (2003). 



M.I. D'yakonov, and V.I. Perel, Sov. Phys. Solid State 13, 
3023 (1972). 



